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A MATHEMATICAL MUNCHAUSEN 
Paul Schillo 


At the time of this writing, ten years have elapsed since the un- 
timely death of I. Neustadt. During this decade, the reform movement 
ii mathematics, which he pioneered, has been marking time — almost as 
if it were afraid to proceed without a leader of his imposing stature. 
But, in order to imbue mathematics with Neustadtian methodology, what 
is needed most is not a leader with Neustadt’s revolutionary zeal but 
a few more revolutionary discoveries like he was able to make. 

A sufficient reason for the failure to discover new mathematical 
truths in the manner outlined by Neustadt is the lack of disciples in 
his school. Every reformation must expect the forces of reaction to 
have their innings, but, in this case, the game is being forfeited 
without contest. Most mathematicians have an erroneous conception of 
Neustadt’s position, some have had no occasion to use his famous | aw, 
with the result that they give him little thought, and instances are 
known where otherwise mature mathematical scholars have failed to 
acquaint themselves with even the rudiments of his philosophy. There- 
fore, while this article is simply a brief account of the man and his 
major work, it is hoped that it might serve to correct some of the 
misconceptions of Neustadt’s intent and so help to proselyte his 
program, 

One notion has it that Neustadt sought to destroy mathematics. 
Nothing could be further from the truth; his orientatim@ was positive 
at all times, for what he wanted to see was its reconstruction along 
modern lines. His sole aim was to make the structure of mathematics 
functional, to have it adopt scientific method and, thereafter, take 
its place among the modern sciences. Because of the cruel way in which 
hasty critics have distorted his progressivism, it is a wonder that 
he has not been accused of having suggested the tearing down of the ca- 
thedral of Notre Dame before the erection of the United Nations’ building. 

A good example of the way in which he has been misunderstood is 
the popular interpretaion of his facetious disproof of Malthus’ Law. 
From the readily verifiable fact that two can live as cheaply as one, he 
proceeded to show by conventional mathematical methods that ten billion 
can live as cheaply as one. To naive critics, this might seem to be an 
attempt to discredit mathematics as a whole; actually, it was intended 
merely to show the absurdity of mathematical induction, a technique of 
which he strongly disapproved. He had no quarrel with fundamental, 
time-tested devices like addition and multiplication and used them 


frequently in his experiments. 
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[t was not without justification that Neustadt came to the con- 
clusion that mathematics was suffocating in its steadily worsening 
state of introversion, He felt that it had accepted the egocentric 
predicament centuries ago and had given up trying to extricate itself. 
To him, it was autopophagy, all this waste of cognitive energy in pro- 
ducing theorems, which were no more than reworded statements taken 
from the particular axioms that happened to be in vogue. “The present 
pace of mathematics is that of a snail,” he once said. “And what is a 
snail? It’s just a snake that has swallowed its tail” It was this 
irritation with all self-centered philosophical positions that goaded 
him on to his classic proof of the existence of a third personal re- 
ality,. the essence of which is found in his often-quoted remark: “Er 
iszt; daher ist er!” He was sure that mathematics had been feeding 
off its own fat for so long that it had become easy prey to the pre- 
datory natural sciences which picked at its bones like it was cold 
turkey. He envisioned its rebirth and hoped to teach it how to hunt 
for new truths in order to survive in a fiercely competitive age. 

Unfortunately, the complete story of Isack Neustadt has yet to be 
written. The information about his early life is fragmentary. If he 
wrote anything prior to the formulation of the great law which bears 
his name, there is no record of it. Only a sketchy description of his 
formative years has been gleaned from the close-mouthed villagers of 
Dilidorf, Pennsylvania where he was born and raised. His buddies in 
the Grand Amy of the Republic, in which he served as a private until 
he suffered a head wound at the battle of Gettysburg, have added little 
to the portrait, except to agree that he kept to himself and was af- 
fectionately referred to as “Pop”. 

Concerning his life in Dilldorf after the war, this much is known: 
The little indications of budding genius caused the villagers to con- 
clude that he had been hexed, a consequence of having taken up arms 
and of having been encamped for a long time near the wicked city of 
Philadelphia. He began to read non-religious books, some of them in 
English and some of them concerned with the heinous heresy of science. 
Also, he showed no inclination toward physical labor and did nothing 
but pick a little fruit during the harvest in order to get money for 
books. It is to the lasting credit of the gentle people of Milldorf 
that, despite their ignorance of I sack Neustadt’s greatness, they saw 
to it that he did not want for anything essential to his survival. He 
was given shelter in the smithy and shared the victuals in the biack- 
smith’s dinner pail. In fact, he lived more comfortably than any of 
his benefactors. But the need that these good people failed to satisfy 
was his greatest need, the hunger to be understood. The admirable 
patience, with which he endured those pangs throughout his early years, 
was to be given good use later when he had to face the hard fact that 
most mathematicians were simply Dilldorfers with academic degrees. 

His escape from the subtle tyranny of his idyllic, native village 
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was one of those remarkable little accidents which, taken all together, 
play the dominant role in the pageant of science. During the roaring 
twenties three gangsters out of Cicero, Illinois rubbed out a Wet 
lobbyist in Washington. The federal men soon knew who the hoodlums 
were 48 well as their devious getaway route via the backroads of five 
states. But, until they spoke to Isack Neustadt, none of the witnesses, 
who had seen the black limousine on its chicken-killing, wash-dusting 
journey through the hinterland, could help them refute its owner’ s 
claim that it had not been outside of the Chicago environs since its 
purchase. And it was only by a stroke of luck that they were able to 
contact Neustadt. For, due to their dislike of the government and its 
minions, the people of Dilldorf treated the investigators coldly; and 
those who spoke English maintained with righteous indignation that, 
during the fruit harvest none of them could possibly have been so 
idle as to gape at a “devil’s chariot”. But, just as the government 
men were about to drive off, a village wag made some remark about “der 
verdammte Faulpelz”, which caused the dour Dilldorfers to laugh most 
heartily. Questioned about this outburst, an unusually helpful young 
interpreter explained that they were amused by the thought of one who 
might not have been so busily engaged at the time of the vehicle's 
passing, and that was I sack Neustadt. 

As it happened, Neustadt had a vivid recollection of the incident. 
In keeping with his neighbors’ opinion of him, he had been napping 
under an apple tree alongside the road. He told of having a strange 
dream about the war, a part of his life about which he had all but 
forgotten. A falling apple struck him on his temporal scar and he 
awoke in the same state of intense alertness as when he had regained 
consciousness in that gully at Gettysburg. It was at that moment that 
he saw the car pass by in a swirl of dust. He could not make out any 
of its details, not even its color, but its license plate seemed to be 
suspended motionlessly before him. It was an Illinois plate with the 
number 2-]12-21-2. Even if he had not been endowed with a photographic 
memory for numbers, he would have remembered that particular one, 
because he had a profound respect for numerical palindromes. 

Since his testimony was critical in the successful prosecution 
of the three torpedoes, he received a five thousand dollar reward 
which, after taxes, was reduced to $2,122, 12 by a remarkable coincidence. 
However, by going to court and by accepting money for just remembering 
a number without doing any laborious reckoning, he became an anathema 
to the people of Dilldorf. Therefore, he hiked to Harrisburg and hopped 
a freight train on itsempty run to Florida, He took up residence in 
a boathouse on the ocean side of ‘liami Beach, where, because of his 
frugality, he was able to live out his life comfortably on the re- 
ward money. It is a matter of public record that he never received 
remuneration for his many published works; money meant nothing to him. 

In the stuffy refrigerator car which took him south, he began to 
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contemplate his great law. Meditating on the events which had altered 
his life, he considered the cause and affect relationship that had 
to exist between the falling apple and the remarkable number 2- 12-2 1-2, 
He thought of another falling apple in the past and of the man who 
had seen its significance. After all, was not the latter’s name merely 
the anglicization of his own? A combination of his native mysticism 
and the semi-delirium, induced by the stifling atmosphere of the car, 
led him inevitably to the conclusion that he was the reincarnation of 
Isaac Newton and that there was something of great importance to be 
discovered in the number 2-]12-21-2. Whether or not one is disposed to 
believe that there was anything supernatural about such an accumulation 
of coincidences, it must be agreed that the subsequent behavior of 
Isack Neustadt constitutes an excellent illustration of the James - Lange 
theory. For, even if he were not Isaac Newton, he acted enough like 
his prototype to make any distinction seem irrelevant. 

Both literally and figuratively, it can be said that Neustadt’s 
Law was incubated in the warm sand of the Atlantic Ocean, for, during 
his first year of peripatetic study on the beach, he wrote the number 
]-]2-21-2 over and over again on that golden slate. When the first 
burricane in his experience seemed about to rip the boat-house from 
its pilings, his most brilliant idea struck with greater intensity 
than the wind. He converted the dashes in 2-12-21-2 to minus signs 
and put plus signs before all but the leftmost of the other digits, 
obtaining 2- 1 +2- 2 +1 - 2, a finite series with zero sum. He 
knew that other six-term series had the property of adding to zero, 
but he was not to be dissuaded from the further consideration of this 
particular one. 

Within a month from that fateful moment, he made the final altera- 
tion in his series, whence it became .2 - 1 + 2- 2+ 1 - .2. which 
as the reader knows, is now called“ Neustadt’s basic series” or,by 
some writers, simply “Neustadt’s base’. According to his own account, 
the first arithmetic progression which he used was 1,2,3,4,5,6. When 
he produced the derived series (or derivative) .2- 2+6-8 +5- 1.2 
by multiplying the successive terms in the base by the respective 
terms in the progression and observed that it too had a zero sum, he 
was so overwhelmed by emotion that he had to sit down beside a palm 
tree and there fell into a state of unconsciousness. 

However, with the discovery that, for any six-term arithmetic 
progression, the derivative of his base had zero sum, his work had 
just begun. He was then faced with the herculean task of proving that 
his observation did indeed have the status of a mathematical law. 
This he accomplished in seventeen months without any outside assistance. 

Before beginning the experiment which was eventually to prove his 
law, he took as his control group twenty zero sum bases different 


from .2- 1+ 2- 2 +1- .2 and listed two thousand distinct arith- 
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metic progressions to be used on the bases. Among these progressions 
there were many like 2.71828, 3.14159, 3.56490, 3.98821, 4.41152, 
4. 83483 so it should be evident that he did not limit the generality 
of his results by selecting special progressions which would make 
his work easier. His ability to calculate mentally helped him to 
complete this experiment in a shorter period of time than it would 
have taken ten ordinary mathematicians. But, despite this demonstra- 
tion that an exceptional person can undertake an immense mathematical 
project and, single-handedly, bring it to a successful conclusion, it 
ought to be clear that the search for mathematical truths can best 
be carried on in well-equipped laboratories by highly trained staffs. 
In the future machines and teamwork must be at the heart of mathematics 

or it shall forever remain an agglomerate of results from uncoor- 
dinated individual effort. No doubt, those years spent in solitary 
labor on his famous law led Neustadt to organize the first mathematical 

cooperative and to preach the doctrine of centrally controlled research 

as an alternative to the inherently feeble laissez faire mathematics, 
But his hopes that this small group would be the nucleus of a uni- 
versal organization were based on the faulty assumption that others 
were as devoted to this single objective as himself. The cooperative 
lasted less than four months and did nothing constructive with its 
modest resources except to contribute to the Lemke-for-president 
campaign fund, 

On the first day of his experiment to show that, for any six-term 
arithmetic progression, the derivative of .2-1+2-2+ 1- .2is 
itself zero sum, Neustadt had a fortunate accident. Using the arith- 
metic progression -3,-1,1,3,5,7 on his base, he had just obtained the 
zero sum derivative -.6 + 1+2 -6+5- 14, his 45th verification 
of the law without a single contradiction. He was getting tired and 
mistook this derivative for one of the bases in his control group. 
When he applied the progression -3,-1,1,3,5,7 to -.6+1+2-6+5- 14 
he found that its derivative 1.8 - 1 + 2 - 18 + 25 - 9.8 was zero 
sum, something that had not happened previously to a derivative of a 
base in his control group. And, the next morning, when he realized 
his mistake, he discovered something more amazing. By applying the 
progression 1.6,1,.4,-.2,-.8,-1.4 to 1.8 -1+2 - 18 + 25 - 9.8, a 


derivative of a derivative (or a second derivative) of his base, he 


2.88-1+.8+3.6- D+ 13.72. 


obtained a third derivative with zero sum: 
He immediately revised his program for the experiment in order to 
prove the law as itis formulated today: For any arithmetic progressions, 


all derivatives of any order of the base .2- 1+ 2-2 + 1- .2 are 


zero sum. 

Right after this law was copyrighted on Flag Day in 1932, many 
mathematicians thought that it could be proved much more easily by 
their conventional devices than by Neustadt’s scientific method. How 
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wrong they were is evidenced by the fact that, despite numerous attempts, 
none of them have been able to come up with a deductive proof of it. 
Their sole contribution has been a forceful demonstration of the fact 
that the time Neustadt devoted to his successful project is negligible 
when compared with the time used in their futile efforts to belittle 
his methodology. 

Now, unlike many mathematicians, Neustadt was no absolutist.He 
knew that the real numbers were phenomena like atoms and organic cells 
and that the laws which govern their behavior can be ascertained only 
to a degree. This degree of certainty, or truth value, of a known law 
must be quantified by assigning to it a number between zero and unity, 
commonly called its probability. He determined the probability of the 
truth of his law to be a number which differed from unity well beyond 
the thousandth decimal place, which is so close to certainty that, 
should it ever fail, one might do well to look for some catastrophic 
alteration in the universe. The manner in which he obtained this pro- 
bability from his experimental data depended on a special application 


of the famous formula: 


which Euler communicated to Diderot. As if anticipating the current 
popularity of electronic computers, he used the binary system because 
of its simplicity and, consequently, set a = -2 and 6 = 2, He let v 
be the number of experimental verifications of his law andc the 
number of contradictions. Then, n was made | or the vth prime according 
as v was or was not zero. After he evaluated the integer x, the re- 
quired probability was determined by the empirical formula: 


P ore .. Or cern 


Because his experiment yielded no contradictions of the law, this 


formula was reduced to the simpler formula: 
Pegs sy 


and P approached | asymptotically as the experiment progressed and 
the number of verifications increased. 

With the completion of this work, Neustadt did not rest on his 
laurels. But it is not the intention of this article to go into all 
that he did. It is well known that, as a champion of science, he led 
many crusades. Yet, the power of leadership failed to corrupt him and 
he retained throughout his life that golden streak of humility which 
is the banner of true genius. He was as cordial with the mountain folk 
of central Florida as with the polished aristocrats of suburban New 
York. So unobtrusive was he in his frequent contacts with people along 
the beach that, unless they had gone to Miami with the intention of 
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meeting him, they were unlikely to realize how important a person he 
was. This aura of simplicity, which made him such an amiable man, in- 
fluenced even his most profound thoughts. It enabled him to see that 
complex numbers were just ordinary numbers perceived by complicated 
people. 

It is vain but interesting to contemplate what new truths I sack 
Neustadt might have uncovered had he not met with sudden death in 
his favorite retreat by that seaside coconut palm, under which 
he would sit by the hour like the great Gautama. Who knows what secret 
of the universe was locked up in that last message which he left in 
the sand? Perhaps someday a devoted disciple of his will look at that 
eleven term mathematical legacy: 


.2-4 + 27 - 96+ 210 - 302.4+ B4- B2+ B81 - 20+ 2.2 


and will see its significance. But, as a sagacious humorist put it at 
a recent testimonial dinner, it should be hoped that this favored 
fellow will be lying under an apple tree and not a coconut palm when 
he gets his inspiration; otherwise, he might not live to develop it. 


Buffalo, New York. 





Electrical Analogies 


a new field for mathematicians and physicists 


In applying electrical analogies to various physical phenomena, 
mathematicians and physicists at Lockheed have developed a highly 
effective method of solving complex dynamic and static problems. 


They have been establishing the analogies through use of a Network 
Analyzer so advanced in concept that it is one of, perhaps, three or four 
instruments of its type in the country. Their activities involve solving 
partial differential equations by substituting finite difference equations 


Through these techniques, they have been able to solve problems 
that are, at best, extremely difficult to solve on even the largest 


and most modern digital computers. 


Most of Lockheed’s endeavor in this field lies in unexplored areas. 

The dynamic and static problems are on a scale never before approached 
Both creativity and ability to exercise individual initiative are required 
because of the frontier nature of the work in evolving the analogies and 


relating them to a wide range of phenomena. 


It is a new and ideal field for those possessing a physics major and 
mathematics minor or mathematics major and physics minor. 


The growing importance of this area has created a number of new 
positions. Because this type of work is so recent in origin, experience is 
not required. Inquiries are invited, Address E. W. Des Lauriers 

at Department 90-15. 
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THE ISOPERIMETRIC PROBLEM IN THE PLANE 


Glen E. Bredon 


1, We shall consider here the isoperimetric problem for convex 
regions in the plane and its extension to more general sets. Part | 
is completely elementary as regards the methods used and presupposes 
next to nothing on the part of the reader. In it we prove the existence 
and uniqueness of the solution to the isoperimetric problem in the 
plane, restricted to convex regions. The uniqueness of a solution, if 
one exists, can quite easily be extended to regions bounded by a Jordan 
curve, but the existence becomes a little more difficult to show. In 
the second part we consider the existence proof for these more general 
regions, and also formulate and prove a very general isoperimetric 
theorem valid for arbitrary measurable sets in the plane. This part 
presupposes a very little knowledge of point set topology and measure 
theory, and a few statements are left unproved. Some of the places 
where the words “clearly” or “obviously” are used are admittedly actual 
gaps in the exposition. The author feels, however, that these gaps can 
be easily filled in by the intelligent reader by straightforward, if 


perhaps long, reasoning. 
Part I. 


2. A convex region means here a closed set R in the plane such that 
if p and q are in R, then the line segment joining them lies in H. 

The isoperimetric theorem for convex regions in the plane says that 
among all convex regions in the plane of a given area A, the circular 
disc of area A, and only this disc, has minimum circumference. Clearly, 
this is the same as saying that the circular dise has the largest 
area among convex regions with fixed circumference, and again, these 
are equivalent to saying that for any convex region of area A and cir- 
cumference C, A S C*/47, with equality if and onlyif the region is 
a circular disc. (A and C are understood to be finite. ) 

We note here for completeness that the length of the boundary of 
a convex region, and indeed, in a certain sense, the length of any 
curve, 1s defined as the least upper bound of the lengths of the in- 
scribed polygons. Note that if C is finite then the lengths of all 
the chords, and hence of those from any particular point, are bounded 


and so the convex region is itself bounded 


¥ ‘ 
Lemma 1. If y andy’ are convex curves (i.e. curves that enclose convex 
, by . 
regions) with lengths C and C’, and if Y surrounds y’ then C eC 


Proof. By the definition of C' we may assume that y’ is a convex polygon, 


The proof in this case is effected by successively replacing arcs of 
produced 


y by line segments which are just the sides of the polygon y’ 


until they meet y (see fig. 1). This operation clearly decreases cir- 
the circumference of 


cumference and we finally arrive in tiie end at 


‘ 
y', giving us our result. 
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Figure 1 


This lenma shows, among other things, that every bounded convex 
curve lias a finite length. (That is, it is “rectifiable’.) We need 
only enclose the curve in a suitable circle. 

Next we develop briefly the interesting operation of symmetrization, 
due to Steiner. (See (5) and also Pélya (4).) 


3. Given a bounded convex region & and a line Ll, the “axis of symmetri- 
zation’, we define K, the symmetrized image of Kk with respect to l as 
follows. Consider lines r which are perpendicular to | and cut H, Let 
r' be the segment of r which is the intersection of r with R. Let r, be 
the segment of r which has the same length as r‘ and has its midpoint 
on l, We then define K, as the union of all these line segments r ). 


(See Fig. 2). It is easy to see that R, 18 again convex. 





Figure 2 


Lemma 2. Area (hy) Area (R) and Circ (R,) $ Circ (Rk). 


Proof. That the area is invariant is an elementary application of 
integral calculus. We reduce the second statement to the corresponding 
one for polygons, Let P; be a polygon symmetric to | and inscribed in 
KR, such that Cire (P)) 2 inte (kK) - €, where € is anarbitrary positive 
number, Let P be the polygon inscribed in H whose vertices are mapped 
by the symmetrization into the vertices of P,. P, is then, in fact, 
the symmetrized image of P with respect to l. If the theorem is true 


< 


for P and P, then we have, using Lemma 1, Circ (k,) $ Circ (F)) 4+ ES 
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Circ (P) + € $< Cire @) + € and, since € is arbitrary, the theorem 
would be proved. The theorem for polygons clearly reduces to showing 
that the sum of the lengths of the lateral sides of a trapezoid with 
fixed bases and altitude 1s least when the trapezoid is isosceles, 
which in turn reduces to showing that among triangles with common base 
and given altitude the isosceles one has the least perimeter. This 
latter can be seen by use of an ellipse whose foci are at the ends 
of the base of the triangles and whose minor semi-axis is equal to 
the common altitude of the triangles. Since the sum of the distances 
from any point on an ellipse to the foci is constant we have (see 
fig. 3) AC + CB 2 AD + DB AE + EB giving our result, 


Figure 3 


Lemma 3. If R is a bounded convex region and not a circular disc then 
there exists a centrally symmetric region /?' which is convex, not a 
circular disc, of the same area as Rk and of no greater circumference 
than Rk. 


Proof. We note first that any figure F, all of whose synmetrized images 


are circles, is itself a circle. If suffices to show that, for any n, a 


regular 2” gon can be inscribed in F. Select an arbitrary axis of 


symmetrization l and let P, be the chord of F which goes into the 
diameter of ¢, perpendicular to |. This is our 2 - gon. Say that a 
yne-l 


- gon Pe has been inscribed in Ff. Then, given a side S of ee 


1 
consider the two symmetrizations in lines respectively perpendicular 
and parallel to S. The first shows there is a point of the boundary 
of F at the correct distance from S, and the second shows that this 
point must be precisely in the correct position to be a vertex of a 
.” This discussion is 


regular gon builton the vertices of - 


admittedly vague and the reader is urged to ahi a drawing and con- 
vince himself of the validity of the argument. 

We now return to the region R in the statement of the present Lemma, 
By the above we can assume that R 1s symmetric with respect to some 
axis l, since we can always find an axis of symmetrization which doesn’t 
send K into a circle. Now diviae K into two parts by this axis of 
symmetry and flip one side over as in figure 4. This clearly turns R 
into a centrally symmetric region with the same area and circumference, 


and which 1s convex and not a circular disc. 
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4. Proof of the [soperimetric Theorem. 

Theorem 1, (Uniqueness) If among all convex regions of circumference 
C, h has a maximum area, then R is a circular disc. 

Proof. Using lemma 3, we know that if A is not a circle then there 
must be a convex, centrally symmetric region R’ with the same area and 
circumference as K which is also not a circle, Let g be the greatest 
radius of RK’ and s the smallest. Then by using the construction in- 
dicated in Figure 5 we see that s xz, for, if not, we have obtained 
a region with area equal to (g - s)* sin WY + area (R) and the same 
circumference as R, But the only centrally symmetric region with all 
radii equal is the circular disc, which shows that R must be a circular 


dise after all. 


M 


Figure 4 Figure 5 














Theorem 2, (Existence) The circular disc has a maximum area among all 


convex regions with given circumference. 


Proof. Say that K is a convex region with circumference C and area A 
such that C*/A < 47. We can assume R to be centrally symmetric by 
lemma 3 and hence can apply the operation indicated in Figure 5, which 
clearly gives a region which is again convex and centrally symmetric. 
Let HK, be the region which results from n applications of this operation 
and denote by g, and s, the greatest and smallest radii of Rh, respec- 
tively and by Y, the acute angle between them, Each application of 
the transformation leaves the circumference C fixed and increases the 


? R ~ 
area by (g, s,)° sin ¥Y,. Since Rk, is convex and s, the smallest 


radius, i, must lie completely to one side of a perpendicular to s, at 


: > , . ~ > 2. ‘ » j > . . 
its endpoint. Hence we see that s, 2 4, cos, and thus 


2 ae 
) Sin y f= #8 ) / 8 
n n n n n n 


(g -s 


Since the area cannot increase indefinitely (it is bounded, for in- 


2 . ° 
stance, by (C/2)*), we must have (g, -s,)~ 0. Now consider two circles 
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of radius s, and bn inscribedin and circumscribed about R, respectively. 


Recalling lemma I we have the inequalities 
and ms? <A 
n 


n 


which show that s, and g, approach a quantity r with the property 


that C 2a7r and lim A. A' nr. But since A,2A we have that 
n 


4n = G*/A's C7/A < 47, a contradiction, proving the theorem. 


Part II. 


5. We have proved the isoperimetric theorem for convex regions. The 
uniqueness part of the theorem is easily seen to hold for regions 
bounded by a Jordan curve since any such region which is not convex 
clearly cannot be a solution to the problem. We will indicate two 
methods of extending the existence proof to these regions. Both methods 
are concerned with the so-called convex hull ACA )of a region R. This 
is defined as the set of all points lying in triangles with vertices 
in Rk. When K# is closed, (kK) is convex and is clearly the smallest 
convex set containing R. Obviously the existence of a solution to the 
Lsoperimetric problem for the regions we are considering will follow 
if we can show that the circumference of A(f#) is not larger that that 
of Rk. (It is easy to see that the area is actually greater unless R 
1s already convex.) 

The first method of showing that cire (h(#)) S cire (Kt) lies in 
observing that the boundary of A(R) consists of subarcs of the boundary 
B of & together with straight line segments joining points of B and 
also that these straight line segments replace arcs of B joining the 
same points. Anyone fairly familiar with point set topology should 
not have much trouble proving this statement and the existence proof 
of the isoperimetric problem follows immediately. 

The second method rests on an interesting expression for the length 
of an arbitrary rectifiable curve, which we state as a theorem but 
will not prove, although the proof isnot difficult (see Blaschka (1)). 


Theorem 3. For any rectifiable curve y with length L we have 


(*) Tes n(p ,»w) dp dy (Lebes que) 


where n(p,v¥) 1s the number of intersections (possiblyinfinite) of y 
with a line of polar angle ¥ and polar radius op, 

Granting this we have immediately that cire (A(R)) < cire (f) since 
the corresponding n(,/) satisfy this inequality except possibly at 


the set of measure zero where the first of these 1s infinite. 
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6. By use of formula (*) we can state the isoperimetric problem for 
the plane in a very general form as follows. Given an arbitrary measur- 
able set M in the plane, we define the circumference of M by the 
integral (*) where n(p,¥/) is now the number of intersections of the 
line (p,¥) with the boundary “MM' of ™ (‘= complementation, 
closure). 

Theorem 4, If Mis a measurable set in the plane with area (Lebesque 
measure) A and circumference C, then A < C*/4n, Furthermore, if 
M= int (YW) (i.e. if M is the closure of its own interior) and C is 
finite, then equality holds only if M is a circular disc. 

Proof. Since the inequality is trivial in case C is infinite, we will 
assume C to be finite. Given M we will make a series of replacements 
which do not decrease the isoperimetric quotient A/C? and will finally 
arrive at a convex set by these replacements, from which the first 
part of the theorem will follow by previous results. 

i) Area (> Area (MW) and WN MCMNM' go that Cire (M) < Circ ) 
We let K = MW, 

ii) Let K = int (). (Note that now K = int (K).) This operation 
is easily seen to be the result of dropping a certain set of boundary 
points, so that the circumference is certainly not increased. Say 
that the area is decreased. This implies that there is a set B of 
boundary points for which m(B) > 0 (Lebesque measure). Now for fixed 
b we have m(B) [2° m(p,¥) do, where m(p,¥) is the linear Lebesque 
measure of the set of points in which B intersects the line (/,¥). 
Since this is positive, the set of p for which B intersects the line 
(p,~) in an infinite number of points is of positive measure. Since 
this is true for all ¥, we have that C is infinite, contrary to assump - 
tion. Hence the area of K equals that of R. 

iii) We define a relation ~ on int (K) by: p ~ q (p,q in int (K)) 
if and only if there is no k in K such that p and q lie in different 
components of kK - {k}. This is clearly an equivilence relation on 
int (K) and since the equivilence classes are open sets there can be 
at most a countable number of them, Let {E} be the set of these 
classes and set F E, and F UY F . Note that K ~ F consists of 
boundary points alone and so by the argument under ii) does not con- 
tribute to the area, Now each F, is met by the rest in at most a 
countable number of points and hence Area (K) = Area (F) = = Area (F)) 
and Cire (kK) 2 Cire (F) = > Cire (F,). Let 


a lub {Area (F )/(Cire (F )) ?} 
n n 


(whose existence will follow later). Then we have 


Area (Kk) Area (F)= © Area (F.) <a> (Circ (F))? < 


a (= Cire (F,)) ? a (Cire (F))?2< @ (Cire (K))? 
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with strict inequality in the case that there are more than just one 
F_. Hence there is an Fy such that 


Area (F) Area (F) 
< 


—— §$ 


(Cire (F))2 (Cire (F,))? 


iv) Now consider a single fF. We see that if any line meets F in 
exactly one point, then F, must lie completely to one side of the 
line, and hence by formula (*) the circumference of the convex hull 
n(F,) is less than or equal to that of F,. We have thus arrived at a 
convex set whose 1soperimetric quotient is not less than that of M, 
and hence, by part I, we have proved the first part of our theorem. 
Previous remarks and the essence of the above proof render the second 


part of the the theorem obvious. 


7. I take this opportunity to point out to the interested reader 
some of the different and highly elegant treatments of isoperimetric 
problems existing in the literature. Especially recommended is the 
excellent treatment by Polya (3) and also that of Gourant and Robbins 
(2). Most of the groundwork was done by Steiner (5), who, however, did 
not bother much with the question of the solutions existence. Some 
more advanced material is to be found in Blaschke (1), while Pélya (4) 
gives some highly interesting applications of symmetrization to mathe- 


matical physics. 
REFERENCES 


(1), Blaschke, W., Vorlesungen uber Integralgeometrie, Teubner or Chelsea. 
(2). Courant, R. and Robbins, H., What is Mathematics?, Oxford, Chap. VII. 
(3). Polya, G., Mathematics and Plausible Reasoning, Princeton, Chap. X. 
(4). Polya, G., Torsional Rigidity, Principle Frequency, Electrostatic Capacity 

and Symmetrizatton, Quar. App. Math., Vol. VI, No. 3, (1948), pp.267-277. 
(5). Steiner, J., Gesammelte Werke, Berlin (1881-2). 





The Second in @ Series of Announce 
ments on Progressive Expansion of 
Program and Facilities in Mathematics 
at the Knolls Atomic Power Laboratory 


GENERAL ELECTRIC'S 


KNOLLS ATOMIC POWER LABORATORY 


(SP UitWiNC010H00E) \T\S NOW DOUBLING THE 
STAFF OF MATHEMATICIANS FOR ITS 


MODERN MATHEMATICAL CENTER 


The steadily advancing nuclear program enjoy an atmosphere in which the creative 
at Knolls Atomic Power Laboratory call mind may find its full fruition 
for new and imaginative departures in As members of the Mathematical Analysis 


a F ‘ ‘ om mt ab 
mathematic ranging from the most a Unit, they will participate in the formula 


struse formulations of fundamental prob 
tion of theories to describe new physical 


lems to the digital solution of physical 
situations now being encountered, in eval 


problems. To meet the consequent expan 
uating these theories and adapting them 


sion of its Mathematical Analysis Program 
to numerical solution by digital computers 


the Laboratory plans to increase signif 
ind in evaluating reactor designs. Design 


cantly the number of qualified mathema 
evaluations will focus on the calculated 


ticians now at work here enough new 
behavior of mathematical models and will 


ys ha hee at n fact, t 
openings have n created, in fac ad employ the most modern techniques in 


more than double the present mathema rh 


computer programming nature and 


tical staff, Mathematicians at all degree 
complexity of these operations call for cre 
levels are invited to join this expanding 
atively new approaches and fundamental 
rogran 
prog ’ advances in these techniques. These mathe 
maticians will also have the opportunity t 
As previously announced, a modern build PI , 
deal with basic research in physics, chem 
ing is now under construction, principally 
istry, metallurgy and many other aspects 
for the use of mathematicians and physi R k 
of nuclear science 
cists. This Center will be equipped with the 
finest of facilities, including digital com The program at the Knolls offers the atmos 
puters that rank among the most powerful phere, the equipment, the richness of sub 
available. Here mathematicians, working ject matter and the material benefits con 
both independently and in association witl ducive to a satisfying career in application 


theoretical and experimental physicists, will of mathemati 


A LETTER TO DOR S R ACKER. EXPRESSING YOUR INTEREST 
WILL RECEIVE IMMEDIATE ATTENTION 


Knolls Alomic Power Laboralory 


OPERATED FOR AEC By 


GENERAL @® ELECTRIC 


SCHENECTADY N Y 





STABILITY OF PERODIC 
TIME-VARYING SYSTEMS 


Louis A. Pipes 


Summar y 


This paper presents a mathematical analysis of the stability of a 
class of mechanical systems that contain time-varying parameters 
of a periodic type. The mathematical analysis leads to the solution 
of linear second-order differential equations of the Mathieu-Hlil] 
type. For purposes of illustration the general methods presented are 


applied to the analysis of the motion of the “inverted pendulum . 


I. introduction 


The general theory of oscillatory mechanical systems and their 
electrical analogues whose parameters are linear and constant has 
been extensively developed and is well understood. In recent years, 
considerable attention and effort has been directed to the analysis 
and to studies of the performance of dynamical systems whose para- 
meters vary with the time, Special cases of technical importance in 
mechanical engineering are described in the treatises on mechanica! 
vibrations by Timoshenko [1] and Den Hartog [2]. Time-varying 
systems of importance in electrical Engineering and commnications 
include the induction generator, frequency-modulation circuits, the 
condenser microphone, and the super-regenerator. The theory of these 
circuits is summarized in the treatise by Mclachlan |3). 

The behavior of this entire class of dynamical systems can be de- 
scribed by differential equations that by suitable transformations 
can be reduced to an equation of the following general form, 

d*x 
(1.1) — -Fltx 
d ¢* 


where the function F(t) is a single-valued periodic function of 
the time t of fundamental period T. 

This second-order linear differentia! equation 1s known in the 
mathematical literature as a differentia! equation of the Mathieu- 
Hill type, [3]. The genera! solution of (1.1) depends on the nature 
of the periodic function F(t). If this function can be represented 
by a Fourier series of the form, 


i 8) 


@ 
(3.2) F(ti=A, + ) A cos(nwt ) + ) B sin(nwt ) 


n=1 n=l 
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where w= 27/T, then (1.1) is known as Hill’s equation. If, on the 
other hand, F(t) can be represented by the first two terms of the 
Fourier series (1.2) so that it takes the form, 


(1.13) Pits =a, 


» + A, cos(wt ) 


then (1.1) 1s known as Mathieu's di®%erential equation. 


The classical! solution of Mathieu’s diflerential equation leads 
to linear combinations of Mathieu functions |3!. The classical 
solution of Hill’s equation involves computations with infinite 
determinants [3]. In general, the classical solutions are difficult 
to apply to special problems of practical importance. 

lt is the purpose of this paper to describe two methods for the 
solution of equations of the type (1.1). Une method involves matrix 
miltiplication and the second method is an approximate one originally 
given by lord Rayleigh, [4], but which has been used very little 
in the technical literature of vibrations. Although these methods 
are general, in order to provide a concrete physical example of 
their use, they will be applied to the study of the motion and 
stability of the “inverted pendulum) (2,5). 


II. Equations of Motion of the Inverted Pendulum. 


Consider the inverted pendulum 
of Figure 1}. This pendulum con- 
sists of a weightless rod _ of 
length L to which is attached a 
mass m. The base of the pendulum 
is moved with a vertical accelera- 
tion S by a crank or cam mechanism. 
It will be assumed that the pen- 
dulum undergoes a_ retarding 
frictional force at its pivot point 
0 of a viscous friction type. If 
the angle 9 that the pendulum makes 
with the vertical 1s small, 1t can 
be easily shown 2 that the eq - 
uation of motion of the pendulum is 


- gg - 
(2.2) 6 +2b6- r (g+S)6 = 0 


The term 260 represents the retar- 
ding effect of the viscous frictional 
force of the prvot and g is_ the 
acceleration of gravity. In order 
to remove the first derivative term 


Inverted Pendulum 


Figure 1. 
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5 introduce the following transformation, 
1 ] . 

(2.3) ; 5 - [ot +40?) = 0 
dt” L 


This differential equation is of the Mathieu-Hill type (1.1) with, 
1 is 
(2.4) F(t) = [- (y +S) + 0? | 


Two cases of the motion of the base of the pendulum will now 
be considered. In the first case, it will be assumed that by means 
of certain cam mechanism the base of the pendulum undergoes the 


Ls 


following square-wave acceleration, 
(2.5) S = K, 0<t < T/2 
S =-K,. T/2<t <7 


where K is a constant. 
In the second case it will be assumed that through the action 
of acrank mechanism, the base of the pendulum is given the following 


harmonic motion. 


(2.4) S = a cos(pt) 


The effect of these two distinct types of motion of the base of 
the pendulum on its behavior will now be analyzed in detail. 


III. Motion of the Pendulum Produced by a Square-Wave 
Acceleration of its Base. 


In order to analyze this case, let, 


l l ) 
(3.1) . = (y + K) + b%, k? (gy - K) + b? 


The differential equations of the motion of the pendulum (2.3) 
] I 


now take the form, 
(3.2) Xx = 0< t <T/2 
(3.3) ok weer s 


Let the initial angular displacement and angular velocity of 
the pendulum at t = 0 be Gy and wp respectively. These initial 
conditions in 9 lead to the following initia! valuesiin x(t) by 
means of the transformation (2.2). 


(3.4) x(0) = G5 x, (0) 
(3.5) x(0) (wo+b6,)= x,(0) 
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In order to determine the subsequent motion of the pendulum 
starting with the initial conditions (3.4) and (3.5) it 1s necessary 
to cascade the solutions of di%erential equations of the form (3.2) 
and (3.3) over intervals in time of duration T/2. This can be done 
by the use of a matrix mltiplication scheme [5]. 

In order to do this, the solution of (3.2) may be written in the 
following matrix form, 


x,(t) cosh (kt) sinh (k t)/k, x ,(0) 
(3.4) 
x(t) k sinh (k jt) cosh (k,t) x (9) 


with the notation %, =2 and Xo = x. 
Let 


(3.7) 6, = k 1/2 


then at the time t 7/2, the solution of (3.2) may be written in 


the form, 


a) x ,(7T/2) . cosh (0,) sinh (6, )/k, x, (0) 


x (T/2) k, sinh (6) cosh (6,) x (0) 


The solution of (3.3) may also be written in a matrix form 
similar to (3.4) with ky replacing ky. The initial conditions of 
the solution of (3.3) are the values x, (7/2) and x,(T/2) given by 
(3.8). Hence the solutions for x, and x, at t = T may be written 
in the convenient matrix form, 

(3,9) 
x,(T) c(6,) S(6,)/k, C(6,) S(6,)/k, x, (0) 


x(T) kS(6,) C(6,) k,S(6,) C(6,) x,(0) 
T 
(3.10) GO =k, 9 C(6) = cosh (@), S(@) = sinh (6) 


If the indicated matrix multiplication is performed, (3,9) may 
be written in the following form, 


x,(T) x ,(0) 
(3.11) 
x,(T) x, (0) 


where, 
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cosh (6,) cosh (6,) t (k /k,) sinh (6,) sinh (@,) 
sinh (6,) cosh (6,)/k, + sinh (6,) cosh (6, )/k, 
k ,cosh (6,) sinh (@,) + k, sinh (6,) cosh (O,) 

D cosh (6,) cosh (6,) + (k o/k,) sinh (6,) sinh (6,) 


.The solution at the time t nT or after n complete cycles of 
the variation of the acceleration of the base, can be obtained by 
casading the solution (3.11) n times, It is, 


x, (nT) A Bi” x, (0) 


(3.13) 
x(nT) C dD x,(0) 


As a consequence of the transofrmation (2.2), the angle of the 


pendulum at the time t =nT is given by, 


(3.14) 6 (nT) e7 baT x, (nT) 


The stability of the pendulum therefore depends on the behavior 
of x , (nT) for large values of n. By means of Sylvester's theorem 


in matrix algebra [4], it is possible to write, 


A i 1 (S.,,-DS,) (BS) 


(3.15) a 
D 51 (CS, ) (S.,, 74S.) 


n 


where 


(3,16) S, sinh (an) 


(3.17) 
(A+D) 


cosh(a ——— cosh(6,) cosh{ 0, ) + [k4 + k2/2k k ,] sinh(@,)sinh( 6, ) 


Therefore as a consequence of (3.12) and (3.14) the angle of 


the pendulum at t = nT may be written in the following form, 


e 7 onT 
(3,18) 6 (nT) = —— [x (0S, |, - DS, ) + x, (0)BS, ] 


"s 


Therefore it can be seen that for large values of n the angle 
6 will increase indefinitely and the pendulum will become unstable 
in the vertical position if, 
(A+ D) 

2 


(3.19) a > OT, where a = cosh"! 
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If, on the other hand, 

(3.20) bT 

the angle of the Pendulum 6 will eventually approach zero and the 
pendulum will be stable in the vertical position. In the absence of 
damping, b = 0 and the vertical position will be an unstable one 
unless a is a pure imaginary number in which case the hyperbolic 
functions in (3.18) become trigonometric functions. In order for 
this to be possible, it is necessary that, 


A+t+D 


(3.21) ~— 


If (3.21) is satisfied, the pendulum will perform oscillations 
about the vertical position. Since k, is always positive, it is 
necessary that k, assume negative values in order for the condition 
(3.21) to be fulfilled. This means that the acceleration K must 


exceed that of gravity. 


IV. Harmonic Motion of the Base of the Pendulum, 


A study of the motion and stability of the inverted pendulum 
when its base is forced to execute a harmonic motion of the form, 


(4,1) S a cos (pt) 


by the action of a crank mechanism will be undertaken in this sec- 
tion. The method of analysis is an approximate one suggested by 
Lord Rayleigh, [4]. 

If (4.1) is substituted into (2.1), it is seen that the equation 
of motion of the pendulum in this case takes the following form, 


(4.2) 6+2b6-- [y - ap cos*(pt)]@ = 0 


l 
L 


To simplify the form of this equation, introduce the following 


change in variables, 

(4.3) O(t) = en ¥ty(t) 
and 

(4.4) z = pt 


In the new variables x and y, (4.2) takes the following form, 


1*y 


(4.5) 2 ~ Le/L b*)/p? -a cos (x)/L]y = 0 
x 


Now let, 
(4,4) *=(g/L + b*)/p 


(4.7) > = a/f2L 


( 
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Then (4.5) may be written in the convenient form, 


: d*y 
dx? 


(4.8) [w? ~ 2 cos (x)] y 


This is a form of Mathieu's dif%erential equation and its general 


solution has the form 3, 


(4.9) . aed (ntk)s 


n=-@ 


where j = (-1)* and the a, quantities are constants. k is a para- 


meter to be determined by requiring (4.9) to be the solution of 
(4.8). If (4.9) is substituted into (4.8) and like powers of e/™ 
are equated, the following set of homogeneous equations in the a's 


are obtained, 


(4.10) [(n + k)? 4 w* Ja, - cla, . 


where n runs through all values from — infinity to + infinity, and 


use has been made of the result, 
(4.11) 2 cos (x) ge?” + @°?* 


In order for the homogenous set of equations (4.10) to have 
other than the trivial solution a, 0, for all values of n, it is 
necessary for the determinant of the coefficients of (4.19) to vanish. 
The vanishing of this infinite determinant determines the value of 
the parameter k. The recursion formula (4.10) may then be used to 
determine the a’s in terms of two arbitrary constants. This is an 
outline of the classical Floquet theory of the solution of Mathieu's 
equation [3]. 

If, however, the parameters of the pendulum are such that 


(4.12) 
(4.13) 


as is usually the case, then the following approximate method for 


the solution of (4.8) given byLord Rayleigh [4] may be used. 


Instead of the infinite series (4.9) the following three term 
approximation for y is assumed, 


n=+1 


(4.14) y > a el (ntk)x 
, n 


n=-1 


This approximation is justified because as a consequence of (4.12) 
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and (4.13) the values of the a’s of higher indices are very smal] 
and can be neglected. As a consequence of the smallness of w* and c* 


it can also be shown [4] that the parameter k is also of the order, 
(4.15) k® << 
Hence the recursion formula (4.19) gives the three equations, 


-a + 
=-1 


2 2 
17 (k* +w Jay + 
Cay- a, 
where use has been made of (4.12) and (4.15). In order for this 


system of homogeneous equations to have a non-trivial solution, the 
: ae . ; 
determinant of the coefficients of the @’s must vanish. Hence, 


If this determinant is expanded, the following equation for the 


parameter k is obtained 


(4,18) k? = 2c? ~ w? 


If 2c? > w*, then (4.18) gives two real roots for the parameter 


k of the form, 
(4.19) 2 = + (2c? — w?)? 


Fquation (4.14) then gives the following two linearly indepen- 
dent solutions of (4.8), 
Yp C,e4** {[ce“J* +1 +ceJ*] 
(4,20) ‘h 
Yo Cie °) *(ce"J* +1+ce)*] 


where C, and C, are arbitrary constants. Since the equation (4.8) 
is linear, the general solution is the sum of these solutions and 


it can be written in the form, 

(4,21) y = [Cyed** + Cye"I**) [ce“J* + 1 + ced*] 
This can be written in the alternative form, 

(4.22) y =A cos (kx + p) [1] + 2 cos (x)] 


where A and & are two new arbitrary constants. In the absence of 
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damping, 6 = 0 and the angle of the pendulum may be obtained from 
(4.3) and (4.22) in the form, 


(4,23) P(t) =A cos (wot + p) [1 + 7 cos (pt) ] 


where 


(2.24) Wo E - 


The condition of stability is that the parameter k be a real 
number and hence 


(4,25) 
or 


yee a)? 
(4,26) p > = (2gL) 


That is, in order for the pendulum to execute stable oscil lations 
of the form (4.23) about the vertical, the angular frequency of the 
oscillations of its base must satisfy the inequality (4.26). 

If damping is included, then the pendulum will perform damped 
oscillations about the vertical position of equilibrium of the form, 


(4.27) 6(t) = Ae” ** cos (wot + p) [1 + L cos (pt) } 
provided, ; 
(4,28) p> (2gL + L7b?)* 


where 


(4,29) 


V. Conclusion. 


The two methods presented in this paper for the analysis of 
dynamical systems whose equations of motion are of the Mathieu-Hill 
type are quite general and useful in cases where the classical 
theory is dificult to apply. The matrix multiplication method is 
particularly useful when the variable parameters are of the square- 
wave or saw-tooth type. The approximate Rayleigh method is a 
very useful one when it can be applied and seems to have been neg- 


lected in the technical literature. 
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AN ELEMENTARY APPROACH TO THE USE OF THE 
RATE OF CHANGE CONCEPT FOR SOLVING PROBLEMS 


Raphael T. Coffman 


Foreword 


This paper presents a method for deriving expressions for the rate of 
change of certain algebraic and trigonometric functions and for applying 
the rate of change concept to derive area and volume formulas without, 
at any stage of the development, introducing the use of limits, 

From the time of the ancient Greeks all recognized methods of deriving 
such formulas, except for the simplest, require the use of limits in 
some form, ranging from the Greek method of exhaustion to the integral 
calculus, 

The writer has not found, in searching through such writing on 
mathematics as were available to him, any indication that this method, 
or any other not involving the use of limits, has previously been used 
for solving the type of problems considered in this paper. Some state- 
ments were encountered to the effect that for specific problems no such 


method is known. 


I. Introduction 


A line which moves in a direction perpendicular to itself may be 
considered to generate an area by its motion, the rate of productionof 
area is the product of the length of the line and its velocity, Fig. I. 


Fig. 1 


If the line of motion is not perpendicular nor parallel to the line 
the velocity vector will have normal and parallel components, vx and 
vy Fig. 2. It is apparent that moving a line in a direction parallel to 
itself does not produce area, but only changes the position of the 
line. In Fig. 2 the area swept out by the motion of 1 in the direction 
of vector v in any time will be a parallelogram, the area of which 1s 
1 x altitude. The altitude will be proportional to vz, regardless of 
the value of v, for equal time of motion. Hence: 


1. The rate of area generation by the motion of a line ts equal to the 
length of the line and its normal (perpendicular) velocity. If the 
line is changing in length, one end will have a parallel velocity. 


8] 
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[nis velocity does not generate area, hence cannot affect the rate at 
which area is being generated at any specific time. It determines the 
rate at which the rate of area production is changing, assuming uniform 
normal velocity. 

In what follows, all rates of change will be given in terms of one 
uniform velocity, which may be defined, without reference to limits, as 
change per unit time. If a velocity varies and its magnitude can be 
related to some other quantity, the expression showing the relationship 
sufficiently defines the velocity. Thus if vy taux, vy varies with 
x and vx and is specified as fully as is A in A =7r* for any value 
of r regardless of whether r is a constant or a variable. 

In the problem to follow the proposition (1) given above and the 


following are the basis for most of the proofs: 


2. If two quantities, originally zero, vary in such a manner that their 
respective rates of change are always in a constant ratio their 
magnitudes will always be in this ratio. 

[his 1s true for any uniform rates of change, hence must be true 
if the rates of change vary, for it is true for whatever rate of change 
the quantities may have at any instant and they can change only to 
another value for which it is also true, 

In the notation used the rate of change of a quantity will be in- 
dicated-by v, as vx, v@, etc; where coordinates are used vx will be 
assumed to be a uniform rate of change and vy will be expressed in 
terms of vx e.g. vy of Ux. 

[he area bounded by a curve, the x axis and the ordinate at x wll 
be designated as Ax, similarly with respect to the y axis for Ay. 

Volumes may be considered as generated by moving planes and the same 


arguments apply as given for generation of areas by lines. Hence: 


[he rate of volume production by a moving plane is the product of the 


area of the plane and its velocity perpendicular to the plane. 


II. Problems 


A. Consider the area Ax in the first vy 
quadrant, Fig. 3, under the line 


y x. By geometry this area 1s known q 


to be 4x*. This may be found by the 


rate of change method as follows: Ay 


yux, which since y x kB ZX Ax 











xvy, vy vx, hence vAy XUX 


[The total rate of area production 


9 


a. | 
ls XUxX + XUX 2xvux and the area Figure 3 
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produced is xy x*, hence vx? 


Ax : X ux 


Ax + Ay vAx + vAy Qxux 


Ax = 4(Ax + Ay) 4x? 


B. Next consider the area Ax in the 
first quadrant under the curve 
y x*.The total area A is xy 


xx? x? 


vd 





vAx yux x ?vx 


? 
é 


vAy xvUy, vy Ux 2x ux 


a4 


from problem (A), hence, 


Ax 


vAy = BUR £ ws 











{x vAx x*ux 
3x* ux 


Figure 4 


2 


fhe formula has been derived for the area under the parabola y x 
and an expression obtained for the rate of change of x7, namely 3x “ux. 
by following the same procedure as for y x* the area under the curve 
y x? is found to be % x” and the rate of change of x" found to be 
Ax?vx. Thus it can be shown that the rate of change and the area under 
any curve of the form’ y x" are respectively nx"*! vx and 
[1/(n + 1)] x"*! where n is any integer. 

The same procedure may be used to 
show that the same formula applies 
when nis a fraction of the form 
1/a where ais an integer. 


C. The equation y = /x, Figure a 
is the inverse of y x 
Ax (2/3)y¥?. The same result may 


hence 





be obtained as follows: The rate 
of change of x* has been shown to be 
2xux; x may be regarded as (/x) “, 


2/x v/x from which 





Ux v(/x) ? 








Figure 5 
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yux J/x vx 


XUX 


Vi 


~V% 


l/ f 
2% Vx 


Vx vy 


a 
vAx + vAy = “x vx +4 yx vx (3/2) Vx ux 


vAx Ju 0x 2 
vA (342) Vx vx 3 
{x (2/3)A (2/3) xt? 


f : 7 
Both v(Vx) and the area under the curve are of the forms 


nx”? yx and Ax 


' . 96 
ty following the same procedure for y x3/3 y x02 and so on it 
can be shown that the same yeneral equation is true for any fraction of 


the form n l/a_ (a integer). 
To find the rate of change of a product, such as xy, construct a 


rectangle whose sides are x and y, Figure 6. 

















Ne 


Vx 


Figure 6 Figure 7 


Then vA = yux + xvy = vxy. If Ais constant yux + xvy = 0, since a con- 


stant has zero rate of change. 


1). Find the formula for the volume of a pyramid, cone or any pyramidal 
prism, in terms of the area of the base and altitude. 
The rate of change of a cube is vx? 3x°vx (Fig. 7). The side xis 
increasing at a uniform rate. 


If the side of a square in- 
creases at the same uniform rate 
as the cube in Fig. 7 and has a 
velocity vx perpendicular to its 
surface, the rate of volume pro- 
duction is x’vx and the solid 
produced is a pyramid with base x é 
and altitude x. By giving the 


generatins square the desired com- Figure 8 


¢ 
( 
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ponent of velocity parallel to its surface, the shape of the pyramid 
may be made regular or oblique as desired. A rotary motion in the plane 
of the square xx may be imparted, producing a spiral pyramd., 

The rate of change of the pyramid, x‘vx, is 1/3 that of the cube, 
hence when x of the cube and pyramid are equal the volume of the pyramid 
is 1/3 that of the cube or (i/3)x?. vx 
A 


. 


“Lz 











Figure 9 Figure 10 


In Fig. § an area, A, is shown in the generating plane xx. This area 
has a constant ratio to the area of the square, that is A/x* = constant. 


vA Avx 
Hence ~ constant 57. The volumes produced will be in the same 
Ux x ux 


Volume of solid A A. 9 
THD, Ome > . » from which, Volume of solid A 
(1/3 )x x 


) 


Z*° § 


- A*x where A is the areaof the base and x is the altitude of the solid. 


E. Find the formula for the rate of change of the sine, cosine, tangent, 
etc. 

1. To find the rate of change of the sine and cosine construct, as shown 
in Fig. 10, a circle with radius 1, and tangential vector vx which re- 
presents the rate of change of the angle x. AB = sin x and Ob = cos x. 
The vertical component of vx is cos x vx and, since it 1s the velocity 
with which A is moving upward, it is the rate of change of AB. Hence 
v sin x cos xvx, Simlarly, the vector sin xvx the rate of change of 
cos x and since cos x 1s decreasing this change 1s negative. 


v COS x - S1In X UZ. 
2. To find the rate of change of the tangent and secant, draw Figure 11. 
Here since r 1,AB tan x and UA sec xX. 
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yy Se@C X UX 
The vector AC is the rate of change of 
the tangent since this is the velocity 
with which the length of AB is changing. 
AC 


sec x or AC sec X UV; 


vs 


but v; sec x Ux 
AC sec X + sec x UX sec “xux 
2 


v tan x sec “XUXx 


v, 18 the rate of change of sec x (OA) 


tan x or Vv, tan x Vy 


v sec x tan x + sec xX Ux 








sec x + tan xux, 





3. To find the rate of change of the co- 
secant and cotangent: 
In Fig. 2B, 

OA 


] 





CSC X UX 





rate of change of OA 





Figure 12 


cot x sec xX Ux 


Since v, acts in a direction to decrease OA, the rate of change 
is negative, 
v CSC x - sec x cot x ux 


is the rate of change of AB and is negative because AB is decreasing. 


AC 
AC 


csc x or AC (SC X Uy 
v 1 
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, 
v cot x - CSC X Vy - CSC X CSC X Ux - CSC” X UX 


1. Find the rate of change of sin? x. 


In Fig. 13, make OD AB = sin x, then DG: sin*x, 


The point D has two components of velocity a normal, DE = cos x vx 
Since the rate of change of sin x = cos xvux and a tangential, Di 
sin xvx, since the tangential velocity is proportional to the distance 
from the center, 0. The rate of change of sin’x is the vertical com- 
ponent of the resultant of the two velocities. This component is 


KJ = KE + EJ 

KE = DE sinx =cos x sin x vx 
EJ = DH cos x =cos * sin x vx 
KJ = 2 sin x cos x vx 


This result can be obtained in a less complex fashion by considering 
a square whose sides are sin x. The rate of change of the area, which 
: a Ne . 

is sin°’* is 2 sin x cos x vx. 























Figure 13 Figure 14 


F. In the preceding derivations of volume and area formlas the ratio 


yux ; 
ao has been constant, with the result that the desired quantities 


were obtainable in terms of this ratio. When the quantities are not in 
a constant ratio to each other a modified procedure is necessary. 


In Fig. 14 it is desired to find the area under the curve y = f(x) 
between the values of x = 0 and x = a. 

This area is equal to the area of some rectangle whose base is a 
and whose altitude is d, that is Ax = ad. d then is the average height 
of y in the interval under consideration. (Notation: y = average y.) 

vy b b 


' mek y am a 
The average value of vy is given by — = — or vy = — vx. This does 
vx a a 
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not directly assist in finding the required area though it does suggest 
a method of approach: if another equation, y, = f,(x), for which 
vy; f(x)( ux can be found, the average value of vy, will be the 
average height of y. Let fx) c when x =a and f (2) =h when x = 0, 


vy; y ux 


c-h = 
vy; ney Ux y vx 


ae 


a 


c 


a 
f,€a) - f (0) = Ax 


As an illustration of the method, find the area under any portion 
of one arch of the sine curve from x = 0 to x =a, 
The rate of change of cos x is — sin x vx hence f, (2) = =— COS x. 


A ~- eos a ~ (~cos 0) = 1 - cos co 


If the equation of the curve contains more than one term the area 
under the curve will be the sum of the areas undef the curve of each 
term considered separately, that is if y = f, (x) + f,(x) the area 
under this curve will be Axf,(x) + Axf,(x), since, Fig. 15, the rate 
of area production is the same whether the curves are separate or 


combined, 


J | 
“4 


y > 4) 
vx 
Ax, iew 


|__ wx 
Mm x 

















Figure 15 Figure 16 


G. Find the formula for the volume of a sphere. 
In Fig. 16 the sphere of radius r is cut by a plane p perpendicular 
to the x axis and moving with a velocity vx as shown, The rate of 


: - 2 2 2 2 
volume production is 7 y“vux m(r* —- 2°) vx = 7 (r*ux —-x*vz) 
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and the volume of % the sphere m(r*x — (1/3)x3) {,(2) in the in- 
ternal x =0 to x=r, V/2 m(r?-(1/3)r) (2/3) 7r? and V = (4/3)7 r?, 


vr vz 


~, 





S = Surface 
Q.= Volume 


Figure 17 Figure 18 


H. Find the surface of a sphere. 

If the radius of a sphere, Fig. 17, is increasing at a rate vr, 
the rate of change of volume is S ur since the entire surface is moving 
away from the center with a velocity ur. 


S ur =v. ( (4/3) mr?) tor? ur 


S =47r? 
An alternate solution is as follows: 


In Fig. 18 the sphere of radius r and center O is cut by a plane 
moving with a uniform velocity vx, with A taken as the origin, The 
surface considered will be that produced by the circle of radius h, at 
the intersection of the plane and the sphere, as the plane moves to 
the left in the figure. The rate of change of surface, vS = 27h vz, 


; . vz 
where vz is the tangential velocity. csc zZ, or vz *SC Z UX 
Ux 


and h =r sin z. 


* Qah vz 27r sin z* csc z Ux, 
and 


sin z* csc z ie * wS 2m 1r vx, 
and the surface generated when vx has produced a length 2r is 
2nr-2r =4nr? 
[his also indicates that the surface of the portion of the sphere 


between A and the cutting plane at C is that of the convex surface of 
a cylinder of radius r and altitude AB where 


AB=r = Vet Ww 
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The application of the method is that of integral calculus, except 
for the derivation of the formula for the areas under the curve y = x", 
which 1s direct, and yields the equivalent of the integral of the 
function. 

The method of finding the equivalent of the derivatives of the 
trigonometric functions, which is unrelated to the differential 
calculus method was discovered by the writer prior to his contact 
with differential calculus. 

No attempt has been made to present this material in a form suit- 
table for presentation to high sciool students, but rather to explain 
the basic principles involved and to demonstrate the possibilities of 
the method. |Inverse trigonometric functions were not mentioned, since 
they are closely related to the direct functions. Logarithmic and 
exponential functions are not treated. No success has been attained 
in applying the velocity concept to the derivatim of the rate of 
change of these functions. 


NOTE: The integral of the sine and cosine are directly obtainable 
from Figure 1), sin x being the quantity produced by a velocity of 
cos xvx and | - cos x is produced by a velocity sin xvx inthe interval 
zero to x. By an extension of the method the areas under such curves 
as sin’x, sin’x, x sin x, etc., may readily be found. This method 


was developed after the above paper was submitted, 





Richland, Washington 
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MISCELLANEOUS NOTES 
Edited by 
Charles K. Robbins 


Articles intended for this Department should be sent to Charles K. Robbins, 
Department of Mathematics, Purdue University, Lafayette, Indiana. 


GRAPHICAL CONSTRUCTION TO FIND APPROXIMATE LENGTH 
OF THE CIRCUMFERANCE CF A CIRCLE 


Jerzy Cicero-Pienk owski 























Construct ton: 


Draw a circle. Through its center draw two straight lines (xx) and 
(yy) perpendicular to each other. Draw straight line (x,x,) parallel to 
line (xx) and tangent to the circle. Also draw straight line (y,y,) 
parallel to (yy) and tangent to the circle. This will establish three 
points (A), (B), and (J), and the square, (AOIB) whose sides are 
equal to R, the radius of the circle. From point (B) describe an arc 
with its radius equa] to four times (OA), and intersecting line (xx). 
This will locate point (C). From point (C) describe an are with its 
radius equal to twice ((A) and intersecting line (x,x,). Thus we 
shall get point (J). The distance from point (A) to point ()) will 
be approximately half of the circumferance of the circle. If we re- 
peat the same construction on the left hand side of line (yy) we will 
get point (H) and since (HA) = (AD) 


(HD) = circumferance of the circle - approximately. 
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Proof: 


By construction 


(yy) | (xx) and (y¥,) | (x,x,) 
(y,y,) || (vy) and (x,x,) || (xx) 
he | 1° 1 


(BC) = 4(0A) = 4R 
(CD) = 2(0A) = 2R 


If we draw a straight line through point (C) and parallel to line 
(yy) then 


(CE) = (OA) =R- and YX (AEC) = 90° 


Also we get two right angled triangles: (BCE) and (DCE). 
In A ( BCE): 
(BE)? = (BC)?~ (Ck)? 








(BE) = V (BC)? - (cE)? = V (aR)? - (R)? = 
Vier? — rR? =V15R2 = R V5 ; 
In A (DCE ):; 
(DE)? = (pc)? ~ (CE)? 
(DE) = V(DC)2 — (CE)? = V(2R)? — R? 
Van? —R? = V3R2 = V3. 


(BD) = (BE) - WE) = RY15 - RV3 = R(V15 - V3) 
whence 


(AD) = (AB) + (Bv) = (OA) + (BD) =R+R(VI5 - V3) 


RQ. + VI5 - V3) =R (3.14093....) 


Now the circumferance of a circle equal 27rwhere 7 = 3.14159.... 
Hence one half the circumferance equal r+ 3.14159.... 
Hence (AD) equal one half circumferance of a circle - approximately, 


and (HD) equal the circumferance - approximately. 


Ruffalo, New York. 








IBN HAITHAM’ S PROBLEMS 
AND THEIR GEOMETRIC SOLUTIONS 
A.R. Amir-Moez 


Ibn Haitham, (Alazen) (965-1039), has done a great deal in geometry 
and optics. He solved two problemsby constructing the roots of quadratic 
or biquadratic equations. However, he was never satisfied with his 
solutions. 

1, Problem: The circle (O) with center at O and the point P inside 
(0) are given. Find the direction of the ray through P such that after 
two reflections it passes again through P. 

Solution: Suppose PT is the desired ray (Fig. l). Clearly OT is 
the angle bisector of ZP7Q and TR, the tangent to (O) at 7, is the 
angle bisector of ZM7TQ. Let OP = a, OT = r, OQ= x, OR = y. It is 

2 OQ RQ P x Y¥ - 2 


easy to see r“ = xy and OP = RP’ 1.e€., _— . Therefore we get 
a y+ta 


2ax? +r*x - ar* = 0. Clearly x can be constructed geometrically. 


2. The solution of Ibn Haitham’s problem: Suppose PT is the desired 


ray, (Fig. 2). Clearly TQ, the reflected ray, is perpendicular to OP. 
If TR is the tangent to ©) at T, then 7(POQR) is a harmonic pencil. 


On the other hand 7Q is the polar of R with respect to (0). (on- 


sider the line ML, through P and perpendicular to OP. The lines RM and 
RL intersect (0) respectively at C and B. BM and CL intersect at Q. 
PQOR is a harmonic range, so is LKBR. Therefore MK passes through 0. 
Consider the harmonic pencil O(LKBR) which intersects LM at L, M, R, 
and X. Therefore LPXM is a harmonic range and this shows the con- 
struction. 

Ibn Haitham has solved the following problem by use of biquadratic., 


But the geometric construction of it is desired. 

Problem: Given A and B on the plane of the circle (0) with center 
O, find the ray through A such that its reflection on (0) passes 
through B. 





University of Idaho. 








(Nov. -Dec. 


On A Certain Problem In Mechanics 


Murray R. Spiegel 


In a book on mechanics, a problem similar to the following appeared. 


An object having mass of 1 gram moves in a straight line on 
a horizontal frictionless plane with initial velocity of 
4 cm/sec. If it is acted upon only by a resisting force in 
dynes numerically equal to twice the instantaneous velocity 
in cm/sec, how far from the starting place does the object 


come to rest? 


It was intended that the problem be solved as follows. Letting the 
object move on the x axis starting from the origin with initial velocity 
4 cm/sec, it is clear that the differential equation for the instan- 


taneous velocity 1s 


dv 
{ }) _ -2y 
dt 


subject to the conditions 
(2) 4 and 


Since (1) may be written as 


dv dv 
(3) v o2v or 


dx dx 
we have on integration 
(4) v -~ 2x +C 


Applying conditions (2) it is seen that C 4 so that 


(5) v 4 _ 2x 


Clearly then from (5) v = 0 when x = 2 and so the book states that the 
object comes to rest at a distance 2 cm from the starting place. 

The above problem is slightly misleading in that it is but a simple 
matter to show that the object never comes to rest. To show this, start 
with equation (1) and obtain by separation of variables and use of 


conditions (2) the velocity 


(6) v he-?t 


wiich shows that v is never zero in any finite time and that as a con- 
sequence the object never comes to rest. The solution starting from 
equations (3) hides the time element in the problem and gives the false 
impression that the object does come to rest. 


Rensselaer Polytechnic Institute, 


(Will some reader explain this paradox? Fd.) 











A Proof of the Moment-Area Theorem 


J. P. Hoyt 


In finding the deflection of statically determinate beams with 
uniform cross-sections, the following (“moment-area”’) theorem is 


sometimes used: 


When a straight beam is subjected to bending, the distance 
(referred toas tangential deviation) of any point 8 on the elastic 


curve, measured normal to the original position of the beam 


from a tangent drawn to the elastic curve at any other point A, 
is equal in magnitude to the product of 1/FI and the moment of 
the area bounded by the graph of the bending moment equation of 
the beam, the x-axis, and ordinates of the graph of the bending 
moment equation at A and B, with respect to the ordinate through BR. 
(E is the modulus of elasticity of the material of which the 
beam is made and 7 is the moment of inertia of the area of the 
cross-section of the beam with respect to its neutral axis.) 


Figure I 


In Figure I, let arcABbe the graph of y=f(x), (where f(x), 
f'(x), and f"(x) are continuous throughout the closed interval (a,b)) 
from x = a to x = b, Let AK be tangent to AB at A and AF be parallel 
to () . Since tangent ZEAF = f'(a), EF =(b- a) fla), FD=AC=f(a) 
and BD) = f(b). Hence BE = BD-FD-EF or BE = f(b) - f(a) - (b-a)f (a), 


Now although we have found the exact value of the tangential 
deviation, it 1s instructive as an illustration of summation by,in- 
tegration as well as leading to our desired result to find an integral 
whose evaluation will give BE. 

To do this, divide C)inton equal (not necessary but sufficient) 


subdivisions of length4x by the points x,, x5, %4, «+. ¥,.4- Let P 


be the point on AB whose abscissa is x, and be the point on AB 
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whose abscissa is x;4;. Now x;44 = x; +Ax. Let PR be tangent to 


AB at P and QI be tangent to AB at Q. 
By the formula for tangential deviation given in the preceding 


paragraph, 
BR = f(b) - f(x,) - (b - x) f'(x) and 
BT = f(b) - f(x; + Ax) - (b= x, -Ax) f (x; +Ax) 
so that 
F TKR=BR - BT 


f(x, +Ax) - fx,) -Axf' (x, +Ax) + (b-x,)[f' (x, + Ax) - f' (x;)] 


Now 


for any n, 8O that 


Since f(x) and f’(x) are continuous throughout (a,b), we can ex- 
press f(x, + Ax) as fla) + f'(x,) Nx + O0(A x) and f'(x, + A x) as 
f'( “,) { f"(x)) Ax + o(\ x) where o(A x) stands for a function of Ax 
whose ratio to Ax approaches 0 asAx approaches 0. 

Hence upon substituting and simplifying, 
(b- x) f"Cx)Axs 


And ¥ 


So that 


it } 
b 
lim (b - x,) f"(x) Ax f (b —- x) f"(x) dx 
nom 
a 


t= 


by the fundamental theorem of the integral calculus. 
This last integral can be recognized as the moment of the area 


bounded by y = f(x), y =0, x =a, x = 6, with respect to the line 
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x = 6b. Since the second derivative with respect to x of y inthe ea- 
uation of the elastic curve y = y(x) is 1/EI times M, where Mis the 
the bending moment equation of the beam, the moment-area theorem 
follows. 

Apart from the usefulness of this theoremin the strength of materials, 
it merits inclusion in calculus texts (I have never seen it in one) for 
two reasons: one, in deriving it we use the integration process to 
find the limit of a sum whose limit we already know: two, it can also 
be used to find not-so-easily calculuable moments from easily cal- 
culable tangential deviations. Such problems are also useful to an 
understanding of the theorem. One such example would be: Find the 
moment of the area lying in the first quadrant bounded by the lines 
y = sinx, y=, x= 7/2. This moment would be given by the integral 


, 


n/2 
(n/2 -x) sin xdx = -- 1. 
0 
This involves integration by parts. To apply the moment area-theorem, 
find the simplest function, -sin x, whose second derivative with re- 
spect to x is sin x. Draw the tangent to the graph of y =- sin x at 
x = 0. and find the tangential deviation at x = 7/2. This gives the 


required moment. From figure II, the steps are obvious and simple. 


tan ZDOA = -cos 0 = -1 
AD OA 


AD 
~_| C(w/2,-1) 


D(wia.-"/2) 





Figure II 


U.S. Naval Academy. 





CURRENT PAPERS AND BOOKS 


Edited ty 
H. V. Craig 


This department will present comments on papers poostanedy published in the 
MATHEMATICS MAGAZINE, lists of new books, and book reviews. : 

In order that errors may be corrected, results extended, and interesting 
aspects further illuminated, comments on published papers in all departments 
are invited. 

Communications intended for this department should be sent in duplicate 
to H. V. Craig, Department of Applied Mathematics, University of Texas, 
Austin 12, Texas. 


Principles and Techniques of Applied Mathematics. By Bernard 
Friedman. John Wley & Sons, 440 Fourth Ave. New York 16, N.Y. 315 
pages, $8.00. 

A unique combination of abstract theory and practical application 
marks this book. It is the latest addition to Wiley’s Applied Mathe- 
matics Series, edited by I.S. Sokolnikoff. 

Deploring the widening gap between pure and applied mathematics, 
Dr. Friedman states his belief that this gap is an illusion. He con- 
tends that the study of abstract systems helps solve concrete problems, 
and that the studyof specific problems may, in turn, suggest interesting 
generalizations to the pure mathematician. Principles and Techniques 
of Applied Mathematics therefore sets out to show how the powerful 
methods developed by the abstract studies can be used to systematize 
the methods and techniques for solving problems in applied mathematics. 

Two main themes carry out Dr. Friedman’s approach. The first is a 
demonstration of how the abstract theory of linear operators can be 
used to unify the techniques of applied mathematics. Secondly, Dr. 
Friedman develops and explains specific techniques that can be used to 
obtain explicit solutions of partial differential equations. 

Among the novel features appearing in this volume are’ a presen- 
tation of the Schwartz theory of distributions and the S- function 
with applications; a simple proof for the reduction of a matrix to 
Jordan canonical form, with the proof using generalized eigenvectors; a 
simple technique for obtaining the complete eigenfunction expansion 
for a second-order ordinary differential operator; a formal method for 
solving partial differential equations which includes all transform 
methods such as the Laplace and Fourier transforms; and the useof the 
general theory of the linear operators to obtain the Fredholm theory 


98 
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of integral equations. Admitting that his emphasis is on the techni- 
ques, Dr. Friedman maintains that a logical analysis of a problem 
leads automatically to the proper methods for its solution. 


Richard Cook 





Symposium on Monte Carlo Methods. Edited by H. A. Meyer. John 
Wiley & Sons, 440 4th Ave., New York 16, NY. 382 pages, $7.50, 


The first full-lenth treatment of this subject, the new book is 
the work of twenty-two leaders in this promising field, writing about 
their own research and applications. 

The volume is the outgrowth of the symposium conducted by the 
Statistical Laboratory of the University of Florida, and sponsored 
by the Wright Air Development Center. It is the latest addition to the 
Wiley Publications in Statistics, Walter A. Shewhart and ~~ S. S. Wilks 
Fditors. 

“The possibility of going from a physical process directly to a 
numerical solution without benefit of an intermediary differential 
equation,’ Dr. Meyer states, “is such an unorthodox and revolutionary 
approach that the imagination is kindled to reconsideration of un- 
solved problems of old.” 

Originally refined for solving atomic energy problems, Monte Carlo 
methods have begun to spread into other specialties such as market 
analysis, quality control, public-opinion polling, and genetics. 

The papers appearing in this book range from the simplest to the 
highly theoretical and include the following: An Introductory Note by 
A.W. Marshall; Generat.wn of Pseudo Random Numbers by Olga Taussky 
and John Todd; Phase Shifts - Middle Squares — Wave Equations by 
N. Metropolis; A General Theory of Stochastic Estimates of the 
Neumann Series for the Solution of Certain Fredholm Integral Equations 
by G.E. Albert; Neighbor Sets for Random Walks and Difference Equa- 
tions by Theodore S. Motzkin; Monte Carlo Computations by Nancy M. 
Dismuke; Applications of Monte Carlo Methods to Tactical Games by 
S. Ulam; Conditional Monte Carlo for Normal Samples by Hale F. Trotter 
and John W. Tukey; Monte Carlo Techniques in a Camplex Problem about 
Normal Samples by Harvey J. Arnold, Bradley D. Bucher, Hale F. Trotter, 
and John W. Tukey; and An Application of the Monte Carlo Method to a 
Problem in Gamma Ray Diffusion by Martin J. Berger. 

The balance of the papers include: Stochastic Calculations of Gamma 
Ray Diffusion by L.A. Beach and R.B. Theus; Application of Multiple 
Stage Sampling Procedures to Monte Carlo Problems by A. W. Marshall; 
Questionable Usefulness of Variance for Measuring estimate accuracy 
in Monte Carlo Importance Sampling Problems by John E. Walsh; Expert- 
mental Determination of Eigenvalues and Dynamic Influence Cooefficients 
for Complex Structures such as Airplanes by C. W. Vickery; Use of 
Different Monte Carlo Sampling Techniques by Herman Kahn; A Theoretical 
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Comparison of the Efficiencies of Two Classical Methods and a Monte 
Carlo Method for Computing One Component of the Solution of a Set of 
Linear Algebraic Ej;aations by J. H. Curtiss; A Description of the 
Generation and Testing of a Set of Random Normal Deviates by E. J. 
Lytle, Jr.; Machine Sampling from Given Probability Distributions by 
James W. Butler; A Monte Carlo Technique for Obtaining Tests and Con- 
fidence Intervals for Insurance Mortality Rates by John E. Walsh; and 
Experiments and Models for the Monte Carlo Method by Alwin Walther. 

In addition, an extensive bibliography gives abstracts and sum- 
maries of existing literature, many of them formulated in the authors’ 
own words. A considerable portion of this literature was previously 
available only in the form of government technical reports which were 
frequently of limited distribution. The appendix also contains a sec- 
tion on random digits, a section dealing with empirical sampling, and 
a selection of papers making use of stochastic processes, including 
a bibliography prepared by Dr. David G. Kendall. 


Richard Cook 





Irrational Numbers. By Ivan Niven. John Wiley & Sons, New York 
164 pages. $3.00. 


Number Eleven of the Carus Mathematical Monographs, J rrational 
Numbers was published in August by the Mathematical Association of 
America. It will be distributed by John Wiley & Sons. The author is 
Ivan Niven, professor of mathematics at the University of Oregon. 

This book provides an exposition of some central results on irra- 
tional, transcendental, and normal numbers. There is a complete treat- 
ment by elementary methods of the irrationality of the exponential, 
logarithmic, and trigonometric functions with rational arguments. 

The approximation of irrational numbers by rationals, up to such 
results as the best possible approximation of Hurwitz, is also given 
with elementary techniques. The last third of the monograph treats 
normal and transcendental numbers, including the transcendence of 7 
and its generalization in the Lindemann theorem, and the Gel fond- 


Schneider theorem. 
Richard Cook 





Digital Differential Analyzers. By George F. Forbes, 10117 Bartee 
Avenue, Pacoima, California, Paper cover, approximately 200 pages 
including appendices and nine page index. $7.50 

An applications manual dealing with mathematics in the form speci - 
fically related to digital differential analyzers. A knowledge of 
differential and integral calculus is presupposed. Familiarity with 
more advanced mathematics and with a specific DDA is desirable but 


not essential. 
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The common functions are discussed in some detail, as well as servo 
techniques and limiting methods. Iifferentiation, the more advanced 


functions, simultaneous equations, and implicit algebraic functions 


are described in a more general way. Following a section on errors of 
the DDA, the more advanced sections of the text deal with complex 
functions, coordinate system transformations, conformal mapping, 
partial differential equations, and curve analysis. 

Several illustrative problems are discussed. The MDA is treated 
both as a mathematical device and as an engineering and scientific 


tool. 





This year International Business Machines estimates that 
7,500 Mathematicians will be needed to man Computers now on 
order, of which about 1,500 should be Ph. D’s. About 250 


will graduate this year! 


From Aviation Week, Sept. 17, 1956, p. 25 





PROBLEMS AND QUESTIONS 
Edited by 


Robert E. Horton, Los Angeles City Gllege 


Readers of this department are invited to sumbit for solution problems 
believed to be new and subject matter questions that may arise in study, 
in research, or in extra-academic situations. Proposals should be accom- 
panied by solutions, when available, and by any information that will assist 
the editor. Ordinarily, problems in well-known textbooks should not be 


submitted. 
Solutions should be submitted on separate, signed sheets. Figures should 


be drawn in India ink the size desired for reproduction. 
Send all communications for this department to Robert E. Horton, Los 


Angeles City College, 855 N. Vermont Ave., Los Angeles 29, California. 


PROPOSALS 


285. Proposed by Maj. H.S. Subba Kao, New Delhi, India. 

It is well known that given any two areas enclosed by closed curves 
in a plane there exists a straight line which bisects both the areas. 
Is it possible to construct such a line when the two areas are: a) a 


triangle and a parallelogram, b) two triangles? 
286. Proposed by Joseph Andruskiw, Seton Hall University. 


Evaluate ; 
n/2 : boas 
sin x COs x 


I dx 


P l+cos‘x l+sin‘x 


287. Proposed by Chih-yi Wang, University of Minnesota. 


Show that 
cot(77/11) - cot( 27/11) + cot(47/11) - cot(87/11) + cot( 167/11) Y 11 
288. Proposed by Norman Anning, Alhambra, California. 

A particle moves on a Fermat Spiral, r 2 =a’6, under a central 
force directed to the origin. Show that the law of force involves the 
seventh power of the radius in a denominator. 

289. Proposed by J.M. Howell, Los Angeles City College. 

Prove or disprove the following statement: The diagonal of a rec- 
tangular parallelopiped can be any odd integer other than one or 
five, and the edges relatively prime integers. 


290. Proposed by Leon Bankoff, Los Angeles, California. 


The hypotenuse AC of a right triangle ABC is divided into unequal 
segments x and y by the point of contact of the incircle. Given seg- 


102 
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ments x and y: a) construct the triangle ABC and its incircle, b) with- 


out using the fourth proportional, upon x as a base, constructa rec- 
tangle equivalent to triangle ABC. 


291. Proposed by M.S.Klamkin, Polytechnic Institute of Brooklyn 


Determine the minimum of 


Solutions 


Probability of Convexity of Four Points 
264. { March 1956] Proposed by Norman Anning, Alhambra, California 


four points are thrown at random on a plane. What is the probability 
that they will be the vertices of a convex quadrilateral? 


Solution by Huseyin Demir, Kandilli, Bolgesi, Turkey. The probability 
of three points to form a triangle is evidentlyl. Let then ABC be the 
triangle formed by three of the points. For the points A, B, C, I) to 
form a convex polygon it is necessary that the fourth point /) be in 
the region of ABC where we have excircles. These regions of the plane 


correspond to the angles 4, B, ( of the triangle ABC. Hence 


p (A + B + C)/whole plane = 7/27 =% 
Also solved by Richard K. Guy, University of Malaya, Singapore; 
John M, Howell, Los Angeles City Colleve and the proposer 


A Volume of Revolution 


265. [ “arch 1956] Proposed by Stephen Armstrong, Union College 
Schenectady, New York. 


What is the volume of revolution of the ellipse 


revolved about the line x = y? 


I. Solution by Lawrence A. Ringenberg, Eastern Illinots 


State 
College. Assume a > b > 0 


. Let S denote the sector with central 


angle 
7/2 which is bounded on the left by the lines y 


+x and on the right 
by the given ellipse. If (x,y) is a point in S, its distance from the 
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line y = x is (x -y)/ V2. Since the desired volume V is twice the 


volume generated when S revolves about the line y = x, we may compute 
it as twice the double integral of V27 (x -y) over S. Thus 


V =2V2 7 fflx- y) as 


ae ab 
fly) = . Vp2. yy, A=>———. 
Then , Va? + b? 


0 f(y) A fly) 
v=-2Vv2 rf f (x-y) dx dy+2 /2 of f (x-y) dx dy 
“A “y 0 y 


and svi ned 
y =————__ 

3 Va? + b2 
Comment. If the ellipse is rotated through one-half a revolution, the 
volume V, generated may be computed as the double integral of 
¥27 (x -y) taken over that portion of the xy-plane which lies within 
the ellipse and below the line y = x. It may be interesting to note 
that V, could also be computed as follows, Let the lines y = +x divide 
the interior of the ellipse into four quadrants S,, S., S3, Sy, 
numbered counterclockwise starting with S, = S, S of the preceding 
paragraph. Let V;, u = 1,2,3,4, be the volume generated when S; rotates 
one-half revolution about y = x. Then if V is the volume obtained in 
the foregoing solution and 

4 V2 nab? 

v' —— , obtained from V by interchanging a and 6, we have 


3 Va2+ b2 
V+V' 2 V27 ab Va? + b? 
3 








V,= V/A, V,=V,=V‘/4, and Vo = 


3 


II. Solution by C.N, Mills, Augustana College, Sioux Falls, South 
Dakota. There are two equal generating surfaces revolving about the 
axis y x. The area of each surface is 

ab [ ; b | 
— | - 2 arc sin ———— 
: Va? + b? 
The distance of the centroid of each surface from the axis is 
2/2 a? b 
|r - 2 are sin ——— ] . 
3Va2 4 Va? + b2 


b2 
Applying the theorem that the volume of a solid of revolution is equal 
to the area of the generating surface times the distance through which 


the centroid moves, we obtain 
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4 V2 abn 


3 Va? + 62 


Also solved by Huseyin Demir, Kandilli Kolgesit, Turkey; J.M.C, 
Hamilton, Los Angeles City College; Richard Kk. Guy, University of 
Malaya, Singapore; Sister M. Stephanie, Georgian Court College, New 
Jersey and the Proposer. 


Points Inverse in a Circumcircle 
266.(March 1956] Proposed by Huseyin Demir, Zonguldak, Turkey. 


If M and M’ are points inverse to each other with respect to the 
circumcircle of a triangle ABC then prove that: 


Z BMC + ZBM'C = 2ZA 
Z CMA + ZCM'A 22B 


Z AMB + ZAM'B =22ZC 


I. Solution by Richard K. Guy, University of Malaya, Singapore. 
In triangles COM and 'OC angle O is common and as OV + OM’ = OC? we 


OC Om’ 


have — = Oc" Hence the triangles are similar and JOM'C = YOM. 
In the same way JOM'B = ¥ (BM. Adding these to YOWC and 4OMB we 


have YBMC + XYBM'C nm - XCOM + 71 - XYBOM = XBOC Yee S 
Simmlarly we have {CMA +YCM'A = 24B and JAMB = {AM'B = 2 XC. 


II. Solution by Maimouna Edy, Hull, PQ, Canada. Represent points 
A, B, C, M, M’ by complex numbers Zi» Zoe 2agr 2, z' respectively. 
Let parentheses represent cross ratios and the bars the complex con - 
jugate. We then have: 


() 
This says that the homographic transformation which sends z yp ty 2, 


into 1, 0, © respectively, that is the transformation of the given 
circle into the axis of reals, sends z andz’ into two conjugate complex 


points, 
Now equation (1) reads explicitly 


Evidently 
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Therefore 
4 
a “a z ) ~2 7’ t,-2, +2 “ 
z4°-2 2, z,°2, 
This means that for oriented angles, 


(angle MC, MB) + (angle MWC, MB) = 2(angle AC, AB) . 


The oriented angles form an additive group isomorphic with the multi- 
plicative group of the unit circle. In other words, we may take 
arbitrary measures of our angles and add them (mod 27). The other two 


relations are proven similarly. 
Bankoff’s solution also noted the necessity for proper orientation 


of the angles. 
Also solved by Leon Bankoff, Los Angeles, California; J.W. Clawson, 
Collegeville, Pennsylvania; and the pro pser. 


Limit of a Radical Sus 


267. (March 1956] Proposed by Alan Wayne, Cooper Union School of 


Engineering, New York. 


Va (k - 1)" _ Wn 


If R, — = —— prove that Jam R i. oe 


I. Solution by Paul Schillo, University of Buffalo. Let P(c,u,r) 
denote any polynomial in u of degree t, where 0 < t <r and where 
c is the coefficient of u’. If p is a positive integer and b 7 0, then 
lettingQ.(1/m) denote a function f of m and k such that mf(m,k) is 
bounded for all positive integer pairs (m,k), in which m is greater 
than the greatest integral root of P(b,m,r + p) = 0 andk < m,it 
follows (by the fundamental theorem of the integral calculus) that 


1/p 


a 1 /p a 
k i? 1 
Lim > F ee Ad | = Lim 2 > [4 (-) + a) | 
m 


a-~o mM k 


ae | P(b,m,r + p) =1 


: a 1/ 1/ 
flee] ete" 
. b » + p 


In the given problem, k” - (k - 1)” =P(n, ee n - 1) and m"**= 


P(l, a, n+ VD), so Lin A, ‘c—r39 iS «2 A 


It might be added that the above limit holds if the polynomials have 
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complex coefficients, since it is the line integral of an everywhere 
analytic function. 

II. Solution by the proposer. Make use of the (u, v) rectangular 
plane coordinate system in which u = x + iy, v =x - iy, andi ¥-L 
In this (u, v) plane, the distance between, oo two ~~ P (u,v, ) 
and P,(u,, v,) is given by PP, P, = (uy - u,)*(v, - v,)*. Moreover, the 
omnes of the arc of any réuri tebhe curve v = f(u) joining P, and 


"9 
P= ff (dv/du)” d 
1“ 2 au u. 
dk | 


Now consider the curve v = u", with its points 


P, is given by 


Pi(k/m, k”/m™), (k =0, 1, 2, «22, a). 


Clearly for k>0, P,P, =V ik - kk - 2" / Va *l, Then 


1 
lim R, 2 r f (nu"* 1A du 2VnAn+ 1) 
0 


a~@ 


Similar limits of sequences of series can be evaluated by this 
“ : ” 
transformation” method, 

Also solved by Maimouna. Edy, Hull, P.Q. 


University of Minnesota. 


Canada and Chih-yt Wang, 


Coaxal Circles 


268. [March 19956] Proposed by J.W. Clawson, Collegeville, Pennsylvania. 


Three coaxal circles, centers at A, B and P have the common points 
C and PD. Any straight line is drawn through C cutting the circles 
again in the points L, M and N respectively. Prove that the ratio 


IN/NM equals the ratio AP/PB. 


Solution by Chih-yi Wang, University of Minnesota. Let the feet 
of perpendiculars from A, B and P onto the line LCMN be A’, B’ and P’ 
respectively. Then we have 


CP' = CN/2, CB’ = CM/2, A'C = LC/2; 





P'B'—B'P' A'B' AB’ + B'P' —A'P 
PB ~—SsiBP-—sSCA@BB AB+ BP AP 
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whence LN LN/2 A‘C +CP* A'P’' AP 

NM NM/2 CB’ + PC’ P*B’ PB 
Remark: If the line passing through C is the line CD, L, M and N 
coincide and in this exceptional case the ratio LN/MN is of the form 


0/0. 
Also solved by Huseyin Demir, Kandilli, Bolgesi, Turkey; Maimouna 
Edy, Hull, P.Q.,Canada; Richard K. Guy, University of Malaya, Singapére; 


and the proposer. 
The Product of Power Series 


269. March 1956] Proposed by M.S. Klamkin, Polytechnic Institute of 
Brooklyn, © 
ma 


= n n oJ n a 
- + — 4 — T cee 
Find the sum mA 3! 


n=1 1! 


Solution by Maimouna Edy, Hull, P.Q., Canada. 


@M 


+1 


x” ] "i 

-_ nx” « siteeetins 
» (n.+ ])! x (n + J)! 
n=O n= 0 n=0 


@ 


. git ; 
n ; 


n=0 


This result is valid for |x| < L 


Also solved by L. Carlitz, Duke University; V.D, Gokhale, Atlanta 
University, Georgia; John M, Howell, Los Angeles City College; Alice 
A. Huffman, University of California at Riverside ;Lawrence A. Ringenberg, 
Eastern Illinois State College; Paul Schillo, University of Buffalo; 
Chih-yt Wang, University of Minnesota and the proposer. 


Circles In A Crescent 
270. [March 1956] Proposed by Leon Bankoff, Los Angeles, California. 


A maximum circle is inscribed in a crescent formed by a semicircle 
and a quadrant of a circle. Find a general expression for the radii 
of consecutively tangent circles touching the sides of the crescent, 
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the first touching the maximum circle, the second touching the first 
and so on. 

Solution by Huseyin Demir, Kandilli, Bolgesi, Turkey. Let the 
given circles (0), (0') intersect each other at A and B, and let the 
center and radius of the nth circle be denoted by (0), r, respectively, 

We invert the figure with center at A, k* = AB* = 4R? being the 
power. Under the inversion, (0) is inverted into its tangent line 


BH,, and (0‘) into the line BO’, forming an angle of 2a= 45°. The 


circles ( Q,;), inverse of (0;), from a series of tangent circles in- 
scribed in the above angle. Let ({),) touch BHy at H,. Then we may 
easily find that the radius p, = 1,H, of ({11,) is given by 

Po 


P, 


where 
Po = Qf o = BH tan a= 2R tan a. 
Drawing the common tangent AT,7) to the inverse circles (0), 


(Q.) we write 
Tr, = AT, P,/AT, = AT, + AT. 2, /ATL? =k’p /(AQ? — ni?) 
Denoting the projection of {), on AB by K, we have 
M2 + Ke. o2 
(2k + p,)? +BH? - p? 
4k? + AR pP, + (p, cot a)? 


2 
R Pp, tana 


4R? + 4R p. 4 lo cot’ a 


Substituting the value of P, in the above expression we arrive at the 


desired result, namely 
Ki — 
1 + yl. + sin /8)7" +(1 - sin 7/8)*") sec*" 8 cot /8 
Also solved by J.W. Clawson, Collegeville, Pennsylvania and the 
proposer. 
TRICKIES 


A trickie is a problem whose solution depends upon the perception of the 
key word, phrase or idea rather than upon a mathematical routine. Send us 
your favorite trickies. 
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T 24. A certain physical society is planning a ballot for the election 
of three officers. There being 3, 4, and 5, candidates for the three 
offices, respectively. There is a rule in effect (in order to eliminate 
the ordering of the candidates on the ballot as a possible influence 
on the election) that for each office, each candidate must appear in 
each position the same number of times as any other candidate. What 
is the smallest number of different ballots necessary? [Submitted by 


M.S. Klamkin! 


T 25. (Gaven the sequence 1, 1, 2, 3, 5, 8, ... where @; =@;.,*+ @,. 
find an expression for the nth term, a,» hot in terms of a recursion 
formula, Generalize the sequence so that a, anda, may be any two 
numbers. [Submitted by Ben B. Bowen), 


SOLUTIONS 
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QUICKIES 


From time to time this cnpertaes will publish problems which may be 
solved by laborious methods, but which with the proper insight may be dis- 
posed of with dispatch. Readers are urged to submit their favorite problems 


of this type, together with the elegant solution and source, if known, 


Q 183. What are the ratios of the volumes and surfaces of the inner 


“half” to the outer “half” of the torus generated by revolving the 


2+(y- a)? ré aor [Submitted by V.C. Harris! 


, 


x 


Q@ 184. Simplify 


sin x + sin 2x+ sin 3x +... + SiInnz 


cos x + cos 2x+cos X% +... + cos nx 


[Submitted by M.S. Klamk in| 
Q 185. Prove that 


k ” 7 F 


) fj)" where § 
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r=0 


[Submitted by M.S. klamkin} 
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i(T7) < y(* ~ T): 
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